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RELATIVISTIC PERFECT FLUIDS 
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Abstract. This paper deals with the evolution of the Einstein gravitational fields which 
arc coupled to a perfect fluid. We consider the Einstein-Eulcr system in asymptotically 
flat spaccstimes and therefore use the condition that the energy density might vanish or 
tend to zero at infinity, and that the pressure is a fractional power of the energy density. In 
this setting we prove a local in time existence, uniqueness and well posedness of classical 
solutions. The zero order term of our system contains an expression which might not be 
a C°° function and therefore causes an additional technical difficulty. In order to achieve 
our goals we use a certain type of weighted Sobolev space of fractional order. In [4] we 
constructed an initial data set for these of systems in the same type of weighted Sobolev 
spaces. 

We obtain the same lower bound for the regularity as Hughes, Kato and Marsden [13] 
got for the vacuum Einstein equations. However, due to the presence of an equation of 
state with fractional power, the regularity is bounded from above. 



1. Introduction 

This paper deals with the Cauchy problem for the Einstein-Euler system describing a 
relativistic self-gravitating perfect fluid, whose density either has compact support or falls 
off at infinity in an appropriate manner. 

The evolution of the gravitational field is described by the Einstein equations 

(1-1) Ga/i = Ral3 — -j^gapR = ^T^TajJ 

where gap is a semi Riemannian metric having a signature (— , +, +, +), Rap is the Ricci 
curvature tensor, and R is the scalar curvature. Both tensors are functions of the metric 
Qap and its first and second order partial derivatives. The right-hand side of (1.1) consists 
of the energy-momentum tensor Tap-, which in the case of a perfect fluid takes the form 

(1.2) T"^ = (e + p)M"M^+p(7°'', 

where e is the energy density, p is the pressure and is the four-velocity vector. The vector 
is a unit timelike vector, which means that it satisfies the normalization condition 



:i.3) gapU'^uP = -1. 



1991 Mathematics Subject Classification. Primary 35L45, 35Q75 ; Secondary 58J45, 83C05. 
Key words and phrases. Einstein-Euler systems, hyperbolic symmetric systems, energy estimates, 
Makino variable, weighted fractional Sobolev spaces. 

Research supported ORT Braude College's Research Authority. 

1 



2 



U. BRAUER AND L. KARP 




where V denotes the covariant derivative associated to the metric gai3- Equations (1-1) and 
(1.4) are not sufficient to determinate the structure uniquely, a functional relation between 
the pressure p and the energy density e (equation of state) is also needed. We choose an 
equation of state that has been used in astrophysical problems. It is the analogue of the 
well known polytropic equation of state in the non-relativistic theory, given by 



The sound velocity is denoted by 

de ' 

and the range of the energy density e will be restricted so that the causality condition 
(T^ < 1 will hold. 

The unknowns of these equations are the semi Riemannian metric gap, the velocity vector 
u'^ and the energy density e. These are functions of t and where x"" (a = 1, 2, 3) are the 
Cartesian coordinates on M^. The alternative notation = t will also be used and Greek 
indices will take the values 0, 1, 2, 3 in the following. 

In the present paper we prove the well-posedness of the coupled systems (1.1), (1-2), 
(1.4) and (1.5) under the harmonic gauge condition in asymptotically fiat spacetimes. In 
order to achieve this, we need to rewrite the above equations as a hyperbolic system. 

In astrophysical context the density e is expected to have compact support, or tend 
to zero at spatial infinity in an appropriate sense. It is well known that the usual sym- 
metrization of the Euler equations is badly behaved in cases where the density tends to 
zero somewhere. The coefficients of the system degenerate or become unbounded when e 
approaches zero. It was observed by Makino [18] that this difficulty can be to some extend 
circumvented in the case of a non-relativistic fluid by using a new matter variable w in 
place of the mass density. For this reason we introduce the quantity 



and we call it the Makino variable. 

A similar device was used by Gamblin [11] and Bezard [2] for the Euler-Poisson equations, 
and by Rendall [22] and Oliynyk [21] for the Einstein-Euler equations. 

The common method for solving the Cauchy problem for the Einstein equations consists 
usually of the following steps. 

1. Initial data must satisfy the constraint equations, which are intrinsic to the initial 
hypersurface. Therefore, the ffist step is to construct solutions of these constraints. 

2. The second step is to use the harmonic coordinate condition and to solve the evo- 
lution equations with these initial data. 

3. The last step is to prove that the harmonic coordinate condition and the solution 
of the constraints propagate. That means if they held on a initial hypersurface, 
they hold for later times. 



(1.5) 



p = p{e) = Ke^ , 



K,7 e M+, 



1 < 7. 



(1.6) 



w = M(e) = e— 
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The last step was treated in detail, for example in Fisher and Marsden[10]. The idea 
is to work out the condition VqG"^ = 0. Since our energy-momentum satisfy (1.4) their 
result can be immediately generalised for our case, but for the sake of brevity we have 
omitted the details. 

However, the presence of the equation of state (1.5) introduces an additional step: the 
compatibility problem of the initial data for the fluid and the gravitational field (see (2.10)). 
There are three types of initial data for the Einstein-Euler system: 

• The gravitational data is a triple {M, h, Kab), where is a space-like manifold, h 
is a proper Riemannian metric on Ai, and Kab is the second fundamental form on 
M. (extrinsic curvature). The pair {h,Kab) must satisfy the constraint equations 



• The matter variables, consisting of the energy density z and the momentum density 

appear on the right hand side of the constraints (2.9); 

• The initial data for Makino variable w and the velocity vector of the perfect 
fluid. 

The only type of Sobolev spaces which are known to be useful for existence theorems 
for the constraint equations in an asymptotically fiat manifold, are the weighted Sobolev 
spaces -fffc,5; where G N and 5 e M. These spaces were introduced by Nirenberg and 
Walker [20] and Cantor [5], and they are the completion of C^(]R^) -functions under the 
norm 



Due to the presence of the equation of state (1.5) and the Makino variable (1.6), we 

2 

have to estimate Ww-'^ \\k,5. So it is perhaps worth discussing the estimate of the Sobolev's 
norm of in more details for /3 > 1. For simplicity we discuss this in the ordinary Sobolev 
space = H^{R^). The simplest case is when /3 G N, then < C(||u||Loo)||M||j|^fc and 

there is no restriction on k. When f3 ^N, then we obtain the same estimate, provided that 
k < (3. This bound on k was improved by Runst and Sickel [23] to < /3 -|- |. Applying 
this to (3 = and taking into account the Sobolev embedding < C||M||jf^fc for 

A; > |, we get a lower and upper bound for k: 



The only exception is the case when is an integer. 

Note that for certain values of 7, inequalities (1.8) possesses no integer solution. Hence, 
for these values of 7 it is impossible to obtain a solution to the Einstein-Euler system in 
Sobolev spaces of integer order. So in order to be able to solve the coupled system for the 
maximal range of the power 7, and, in addition, to improve the regularity of the solutions, 
we are considering the Cauchy problem in the weighted fractional spaces Hg^s, where s is 
real number (see Definition 2.1). These spaces were introduced by Triebel [25], and they 
generalize Hk^s to a fractional order. 



(2.9); 



(1.7) 
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In [4] the authors constructed initial data for coupled systems (1.1), (1-2) and (1.4) 
with the equations of state (1.5). This includes the solution to the constraint equations 
(2.9), as well as the solution to the compatibility problem between the matter variable 
{z,j"') and the fluid variables (w,m"), (2.10), in the i/s^^-spaces. Here we will establish the 
well-posedness of Einstein-Euler systems in the weighted fractional Sobolev spaces Hg^s- 

The common way to prove well-posedness is to rewrite the evolution equations as a 
symmetric hyperbolic system. So our first step is to use the the Makino variable (1.6) and 
to reduce the Euler equations (1.4) to a uniformly first order symmetric hyperbolic system. 
This result was announced in [3] and here we present a detailed proof of it. Our hyperbolic 
reduction is based on the fluid decomposition; for alternative reductions see [22]. 

It is well-known that the Einstein equations can be written as a system of quasi-linear 
waves equations under the harmonic gauge condition [6, 7, 26]. The proofs of existence 
and uniqueness either using second order techniques [6, 8, 12, 13, 16], or transferring the 
equations to a first order symmetric hyperbolic system. Fischer and Marsden used the first 
order techniques and obtained the well-posedness of the reduced vacuum Einstein equations 
in and for s > | [10]. This result was improved by Hughes, Kato and Marsden [13], 
who obtained {gap^dtQap) ^ H'^^^ x H'^ for s > |. They used second order theory, and 
took advantage of the specific form of the quasi-linear waves equations, namely, that the 
coefficients depend only on the semi- metric gap-, but not on its first order derivatives. 

Our aim is to prove existence and uniqueness of the reduced Einstein-Euler system (1.1), 
(1.2) and (1.4) with the equation of state (1.5). In addition, we would like to achieve the 
same regularity of the metric as in [13]. But since we have here a coupled system which 
one of them is a first order, the second order techniques of Hughes, Kato and Marsden in 
[13], is no longer available for the present problem. 

In asymptotically fiat spacetimes the initial metric (7^/3(0) differs from the Minkowski 
metric by a term which is 0(l/r) at spatial infinity, and this term does not belong to H^. 
It is therefore more appropriate to consider both the constraint and evolution equations 
in the if^.s spaces rather than the unweighted spaces H^. For the vacuum equations the 
second author obtained well-posedness of the reduced Einstein equations with {gap, dtgap) £ 
if<i+i,<5 X Hg^s+i, s > I and 5 > — |, see [15]. But unlike Hughes, Kato and Marsden [13], 
Karp treated the quasi-linear waves equations as a first order symmetric hyperbolic system. 
The first order techniques have the advantage that they enable, in a convenient way, the 
coupling of the gravitational fields equations to other matter models, in particular, to 
perfect fiuids. In the Appendix we explain the main idea of [15] which allows us to obtain 
the regularity index s > | by means of first order hyperbolic systems. 

A crucial step in the proof of existence and uniqueness of any hyperbolic system is to 
establish energy estimates for the linearized system. In order to achieve this we define an 
appropriated inner-product in the Hg^s spaces, which takes into account the coefficients of 
the linearized system (see Section 5 ). A similar inner-product was used in [15], and here 
we rely on these energy estimates. 

Once we have obtained the energy estimates for the linearized system, we use Majda's 
iterative scheme in order to obtain existence and uniqueness of the quasi-linear symmetric 
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hyperbolic system [17]. This procedure uses the fact that solutions to a linear first order 
symmetric hyperbolic system with coefficients and initial data, are also C^. But 
here we encounter a further difficulty, namely, the right hand side of (1.1) contains the 

2 2 

fractional power see (3.16). So even when w G and w >0, w'<-^ might not be a 

2 

C°° function. We solve that problem by using the fact that that e = w^-^ satisfy a certain 
first order linear equation. Gamblin encountered a similar problem for the Euler-Poisson 
equations [11], but he solved it in a some what different way. 

Our results improve the existence theory of solutions locally in time of self gravitating 
relativistic perfect fluids in several aspects. Rendall studied this problem in [22], but he 
assumed time symmetry, which means that the extrinsic curvature of the initial manifold is 
zero, and therefore the Einstein constraint equations are reduced to a single scalar equation. 
In addition, he dealt only with C^-solutions. In his study of the Newtonian limit of perfect 
fluids, Oliynyk obtained existence locally in time in the weighted space of integer order 
Hk^Si for k > A [21]. Both Rendall and Oliynyk assume that the adiabatic exponent of 
(1.5) satisfies the condition that is an integer. 

The paper is organized as follows: In the next section we define the weighted Sobolev 
spaces of fractional order Hs.s and state the main result. Section 3 has two subsections: 
the first one deals with the hyperbolic reduction of the Euler equations (1.4); in the second 
one we spell out the matrices which describe the coupled equations (1.1), (1-2) and (1.4) 
as a hyperbolic system. In Section 4 we present tools and properties of the if<i,<5-spaces 
which are necessary for us. We also define the corresponding product spaces. The energy 
estimates for the hnearized system are considered in Section 5, there we also define the 
appropriate inner-product. In Section 6 we treat the iteration procedure. A part of the 
steps are standard and known, but some of them require special attention due to the specific 
form of the system (3.24) and the product spaces. In this Section we will use the fact that 
the coefficients of the first order derivatives depend only on the semi-metric gap- Finally, 
in Section 7 we proof the main result. In the Appendix we give a heuristic idea explaining 
how that fact that the coefficients of the wave equations depend only on the semi-metric 
Qap enable us to obtain the desired regularity by means of symmetric hyperbolic systems. 



2. The main results 

We obtain the well-posdeness in the weighted Sobolev spaces of fractional order. So we 
first recall their definition. 

Let {ifjj^^Q C C^(R'^) be a sequence of cutoff function such that, il)j{x) > for all 
j > 0, snpp{ipj) C {x : 2^-^ < \x\ < 2^+^}, ipj{x) = 1 on {x : 2^~^ < \x\ < 2^} for 
j = 1, 2, supp('?/;o) C {x : \x\ < 2}, ipoi^) = 1 on {x : |x| < 1} and 

|9>,(x)| <C,2-H^\ 

where the constant Ca does not depend on j. 
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= 2 and denote the Bessel potential spaces by H'^ 

J{i + \er\m\'dt 

Definition 2.1. For s, 5 G M, 

oo 

(2.1) (ll«b.j' = E2^^^'^''ll(^.^)(2^)lll.- 

j=0 

where fe{x) = f{sx) denotes the scahng by a positive number e. The space Hg^s is the set 
of all temperate distributions having a finite norm given by (2.1). 

The if^ ,5-norm of a distribution u in an open set C M'^ is given by 

M\H,,s{n) = inf ||/||/^^ ,,(ffi3) . 

•'In"" 

Definition 2.2. Let be a 3 dimensional smooth connected manifold and let h he a 
metric on Ai such that {Ai,h) is complete. We say that {Ai,h) is asymptotically fiat 
of the class H^ s if /i € Hf^^[A4) and there is a compact set 5* C such that: 

1. There is a finite collection of charts {{Ui, (pi)}fLi which cover A4 \ S; 

2. For each i, ipJ^iUi) = Er^ := {x G : |x| > rj} for some positive rj; 

3. The pull-back {(pi^:h)ab is uniformly equivalent to the Euclidean metric 6ab on 
for each i; 

4. For each i, {(fi^h)ab - Sab G Hs,s{Eri). 

The ,5-norm on the manifold Ai is defined as follows. Let Uq G M. he an open set such 
that S C Uq and Uq <s A4. Let {xo,Xi} be a partition of unity subordinate to {Uo,Ui}, 
then 

JV 

is the norm of the weighted fractional Sobolev space Hg si-M.). For the definition of the 
norm ||xo'w|lifs(n) manifold see e.g. [1]. Note that the norm (2.2) depends on 

the partition of unity, but different partitions of unity result in equivalent norms. In the 
following we will omit the notation JH, that is, we will write \\u\\h^ ^ instead of \\u\\h^ siM)- 
Since the principal symbol of the fields equations (1.1) is characteristic in every direction 
(see e.g. [9]), it is impossible to solve (1.1) in the present form. We study these equations 
under the harmonic gauge condition 

(2.3) = g^^T% = 0, 

where (?"^ is the inverse matrix of gap- Then the fields equations (1.1) are equivalent to 

the reduced Einstein equations 

(2.4) g^'d^d.gap = Hap{g, dg) - IQ-nT^p + ^Tig^^T^.g^p, 



We restrict ourselves to the case p 
with the norm given by 



\u\ 



where u is the Fourier transform of u. 
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where Hai3{g, dg) is a quadratic function of the semi-metric g^is and its first order deriva- 
tives. 

Since g^'^ has a Lorentzian signature, (2.4) is a system of quasi-hnear wave equations. 
Taking into account the equation of state (1.5), the normahzation condition (1.2), and the 
Makino variable (1.6), then the wave equations (2.4) become 

(2.5) g'^'d^d.g^i, = H^^ig, dg) - 8nw^^ ((1 - Kw^)gap + 2(1 + Kw^)u^up) . 

So the unknowns of the system (2.5) coupled with the Euler equations (1.4) are the 
semi-metric gap, the velocity vector u'^ and the Makino variable w. Note that even if w is 

2 

a smooth function, wt-i might not be smooth in certain regions. 

The initial data consist of the triple {Ai, h, Kab), where is a space-like manifold, h 
is a proper Riemannian metric on M and Kab is its second fundamental form (extrinsic 
curvature). The semi-metric gai3 takes the following data on M: 

goo\M = — 1> gOa\M = 0, gab\M = Kb ^ /, _ 1 o Q 
la I "'1 — i, o. 

— 2'^0gab\M — J^ab, 

The rest of the data ^of^ooi^ are determined by the harmonic gauge condition (2.3). In 
addition, the initial data of the velocity vector u"" and the Makino variable w are given on 
Ai. We denote the Minikowski metric by rjap- 

Theorem 2.3 (Main result). Let | < s < + | and 5 > — |. Assume Ai is asymp- 
totically flat of class Hs+i,s, Kab G Hg^s+i, and (u° — l^u^-^w) |^ G Hs+i,s+2- Then there 
exists a positive T, a unique semi-metric gap, a unite timelike vector u°' and w satisfying 
the reduced Einstein equations (2.5) and the Euler equations (1-4) such that 

(2.7) igapit)-r]ap) e C{[0,T],H,+,^s)nC\[0,T],H,^s+i) 
and 

(2.8) 1,m",u;) e C([0,T],i7,+i,,+2) nCi([0,T],i7,,5+3). 

Remark 2.4 (On the differentiability). Note that we have a lower and an upper bound of 
the differentiability index s, however, in case is an integer, then there is no upper 
bound. 

A necessary and sufficient condition for the equivalence between the reduced Einstein 
equations (2.5) and the field equations (1.1) is that the geometric data {h, Kab) satisfy the 
constraint equations 

..Q^ / R{h)-KabK'^' + {h-'Kab)' = IQttz 

^ ' \ ^^^VbK''^ -^^^V\h^'Kbc) = -^Tif. 

Here R{h) = h"'^Rab is the scalar curvature with respect to the metric h. The right hand- 
side {z,j°') contains the energy density and the momentum density respectively. Thus 
solving the constraint equations (2.9) ensures that the solution of (2.5) satisfies the original 
system (1.1). 
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However, before we solve the constraints, we need to treat the compatibihty problem 
between the matter variables {z,j°') and the initial data for the velocity and the Makino 
variable w. 

This problem can be described as follows: Let denote the projection of the velocity 
vector u°' on the initial manifold Ai and the timelike unite normal vector to Ai. The 
energy density z is the double projection of Ta/s on n° and the momentum density is 
once the projection of T^^/j on n° and once on Ai. Applying these projections to the perfect 
fluid (1.2) results in 

(1 + (1 + Kw'^)) habU'^u'' 



(2.10) 



j'' = w-^ (1 + Kw"^) u''a/1 + habU'^vP 



So the compatibility problem consists of solving (2.10) for u°', when z and j° are 
given. This problem was solved in the iJs^^-spaces and for the same ranges of s and 6 in 
[4, Theorem 2.5]. 

Since, (2;, j") G Hs,s+2, it follows from Proposition 4.10 below that {w,v? — 1,m'^)|^ G 
This explains the difference between the weights for the metric and for the fluid 
variables. We conclude using [1] that there is an initial data set {h,Kab) and (w,m°) in 
the i^s^-spaces which satisfy the constraints (2.9) together with the compatibility problem 
(2.10).' 

Corollary 2.5. Under the assumptions of Theorem 2.3 and in addition under the assump- 
tion that the initial data {h,Kab) and {w,u°') satisfy the constraint equations (2.9) and 
compatibility problem (2.10), there exists a positive T, a semi-metric g, a unite timelike 
vector and w satisfying the Einstein (1.1) and the Euler equations (I.4) for t G [0,T]. 
The regularity of g, u"' and w are the same as in Theorem 2.3. 

3. Symmetric Hyperbolic Systems 

The main result is proved by transforming the coupled system (2.5) and (1.4) into a 
symmetric hyperbolic system. We therefore recall its definition. 

Definition 3.1 (Symmetric hyperbolic system). A first order quasi-linear k x k system is 
symmetric hyperbolic system in a region G gM.^, if it is of the form 

(3.1) L[U] = A''{U)d^U + B{U) = 0, 

where the matrices A°'{U) are symmetric and for every arbitrary U E G, there exists a 
covector ^ such that 

(3.2) io^A^iU) 

is positive definite. The covectors for which (3.2) is positive definite, are called spacelike 
with respect to equation (3.1). 

If ^ can be chosen to be the vector (1, 0, 0, 0), then condition (3.2) implies that the matrix 
A^iJJ) is a positive definite matrix, and we may write system (3.1) in the form 

(3.3) A'^{U)dtU = A''{U)daU + B{U). 
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3.1. The Euler equations written as a symmetric hyperbolic system. It is not 

obvious that the Euler equations written in the conservative form V^T""^ = are symmetric 
hyperbolic. In fact these equations have to be transformed in order to be expressed in a 
symmetric hyperbolic form. Rendall presented such a transformation of these equations 
in [22], however, its geometrical meaning is not entirely clear and it might be difficult to 
generalize it to the non time symmetric case. Hence we will present a different hyperbolic 
reduction of the Euler equations and discuss it in some details, for we have not seen it 
anywhere in the literature. The basic idea is to perform the standard fluid decomposition 
and then to modify the equation by adding, in an appropriate manner, the normalization 
condition (1.3) which will be considered as a constraint equation. 

The fluid decomposition method consists of the projection of equation V JT^^ = onto 
u° which leads to upV^T'^^ = 0, and the projection of these equations on rest pace O 
orthogonal to of a fluid which leads to PapVuT^^ = 0, where Pajs = gai3 + MqM/?- 

Inserting this decomposition into (1.2) results in a system of the following form: 

(3.4a) u''V,e + (e + p)V,u" = 0; 

(3.4b) {e + p)P^pu''V,u(' + PWuP = 0. 

Note that we have beside the evolution equations (3.4a) and (3.4b) the following con- 
straint equation: gapu'^u^ = —1. We will show in Subsection 3.1.1 that this constraint 
equation is conserved under the evolution equation. 

In order to obtain a symmetric hyperbolic system we have to modify it in the fol- 
lowing way. The normalization condition (1.3) gives that Upu'^'V^u^ = 0, so we add 
(e + p)ui3u'''VuU^ = to equation (3.4a) and UaU^u'^'W yU^^ = to (3.4b), which results 
in 

(3.5a) nl'V + {e + p)P%V yul^ = 

2 

(3.5b) To.pu''VuU^ + -^—PWue = 0, 

(e + p) 

where a := is the speed of sound and 

^ap = Pap + UaUp = Qa/i + 2m„M^ 

is a reflection with respect to the rest subspace O. As mentioned above, we will introduce a 
new matter variable which is given by (1.6). The idea which is behind this is the following: 
The system (3.5a) and (3.5b) is almost of symmetric hyperbolic form, it would be symmetric 
if we multiply the system by appropriate factors, for example, (3.5a) by ^ = o"^ and (3.5b) 
by (e + p)- However, doing so we will be faced with a system in which the coefficients will 
either tend to zero or to infinity, as e — )■ 0. Hence, it is impossible to represent this system 
in a non-degenerate form using these multiplications. 

The central point is now to introduce a new variable w = M(e) which will regularize 

2dM 



the equations even for e = 0. We do this by multiplying equation (3.5a) by k^M' 
This results in the following system which we have written in a matrix form: 



K 
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K\e + p)M'P''p 


(3.6) 


cr^ pi/ 

v 





w 



J 

In order to obtain symmetry we have to demand that 
(3.7) M' = . . , 

where k ^ has been introduced in order to simphfy the expression for w. If we choose 
then (3.7) holds and consequently the system (3.6) is transferred to the 



symmetric system 









(3.8) 


\ 





w 

/3 



U 



J 



The covariant derivative takes in local coordinates the form Vy = + T{g'^^ , dgap) 
which expresses the fact that equations (3.8) is coupled to equations (1.1) for the gravita- 
tional field Qaj^. In addition, from the definition of the Makino variable (1.6), we see that 

and a = y/'-fKw. Thus the fractional power of the equation 



W , SO K 



7-1 l+Kw'^ 

of state (1.5) does not appear in the coefficients of the system (3.8), and these coefficients 
are C°° functions of the scalar w, the four vector and the gravitational field gap- 

Now we want to show that of our system (3.8) is indeed positive definite. In order 
to do it we analyze the principal symbol of this system. For each ^T*V the principal 
symbol is a linear map from M x £"3. to M x F^., where is a fiber in T^V and is a 
fiber in the cotangent space T*V. Since in local coordinates = dy + T{g'^^ jdgap), the 
principal symbol of system (3.8) is 









anP'^fs^y 


\ 


(3.9) 




aKP\iy 


(^"e.)ra/3 








\ 




/ 



and the characteristics are the set of covectors ^ for which {C,uA'^) is not an isomorphism. 
Hence the characteristics are the zeros of Q(0 •= '^^^i^.uA'^). 

The geometric advantages of the fiuid decomposition are the following. The operators 
in the blocks of the matrix (3.9) are the projection on the rest hyperplane O, P'^a, and 
the refiection with respect to the same hyperplane, r^^. Therefore, the following relations 
hold: 



7/3 



Or , 



■pa7p u 



po 



and Pp'^Pa = 



/3) 
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which yields 

\ 
/ 

It is now fairly easy to calculate the determinate of the right hand side of (3.10) and we 
have 





/ 1 


\ 








(3.10) 




V 






v 













det 


aKP'^'^i, 








\ 


) 





Since Pg is a projection, 



and since is a reflection, 



det 











pQ7 



det (^?"^r^) = - (det {g^p))-\ 



Consequently, 

(3.11) g(0 := det(e.A^) = -K^det{g^p){u^i,f {{u'^i^f - a^PV^^^} 
and therefore the characteristic covectors are given by two simple equations: 

(3.12) i.u" = 0; 

(3.13) i^^un'-a'P'^pU^ = 0. 

Remark 3.2. The characteristics conormal cone is a union of two hypersurfaces in T*V. 
One of these hypersurfaces is given by the condition (3.12) and it is a three dimensional 
hyperplane O with the normal u°'. The other hypersurface is given by the condition (3.13) 
and forms a three dimensional cone, the so called, sound cone. 

Let us now consider the timelike vector u,y and insert the covector —Uiy into the principal 
symbol (3.9), then 



\ 








\ 





To/? 








/ 



is positive deflnite. Indeed, Fq,^ is a reflection with respect to a hyperplane having a timelike 
normal. Hence, —u^, is a spacelike covector with respect to the hydrodynamical equations 
(3.8). Herewith, we have showed relatively elegant and elementary that the relativistic 
hydrodynamical equations are symmetric-hyperbolic. 
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We want now to show that the covector t„ = (1, 0, 0, 0) is also spacehke with respect to 
the system (3.8). Since P^/jWa = 0, the covector —u^ belongs to the sound cone 

(3.14) i^^un' - ^'P%U^ > 0. 

Inserting t^, = (1, 0, 0, 0) the right hand side of (3.14) yields 

(3.15) (n0)2(l - a^) - aV°- 

Since the sound velocity is always less than the light speed, that is cx^ = ^ < = 1, we 
conclude from (3.15) that t^, also belongs to the sound cone (3.14). Hence, the vector —u^ 
can be continuously deformed to t^, while condition (3.14) holds along the deformation 
path. Consequently, the determinant of (3.11) remains positive under this process and 
hence ti,A'^ = is also positive definite. Thus we have proved. 

Theorem 3.3. Let e be non-negative density function, then the Euler system (1-4) coupled 
with the equation of state (1-5) can he written as a symmetric hyperbolic system of the form 
(3.3), and where is a positive definite. 

3.1.1. Conservation of the unit length of the fluid. 

Proposition 3.4. The constraint condition Qapu'^u^ = —1 is conserved along stream lines 

Proof. Let k{t) be a curve such that k'{t) = u°' and set Z{t) = {u o k)i3{u o k)^, then we 
need to establish 

-Z{t) = 2upVk'[i)U^ = 2u'upV,ul' = 0. 

Multiplying the four last rows of the Euler system (3.8) by u"' and recalling that P^^ is 
the projection on the rest space O orthogonal to we have 

= u'' {Tafsu'^V^u^ + naPW^w) 

= u'^Papu^'Vuu'^ - n!'upVuu'^ + nau'^PWuVO 

= —u'^UpVuU^- 

Therefore, if gapu'^u^ = — 1 on the initial manifold, then it holds along the stream lines 
m". □ 

3.2. The coupled hyperbolic system. In this section we will transform the coupled 
system (2.5) and (3.8) into a symmetric hyperbolic system. We will pay attention to the 
fact that the system will be in a form in which we can apply the energy estimates of [15]. 
That allows us to obtain the same regularity for the gravitational fields as Hughes, Kato 
and Marsden [13] got for the Einstein vacuum equations. Note that our system is slightly 
different from the symmetric hyperbolic system obtained by Fisher and Marsden [10], since 
our system contains a constant matrix as given by (3.21). 

We consider a spacetime (V, gap) of the type M x A^, where is a Riemannian manifold, 
and we denote local coordinates by (t,x''). Set 
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then the reduced Einstein equations (2.5) takes the form 
—g^^dthapo = 

(3.16) 
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{2g'^^daKpo + g^^'daKpb + H^p{g, dg) 

((1 - Kw'^)gap + 2(1 + Kw'^)UaUp)^ 

g'^'dth^fia = g'^'dah^^ 



J/30 

In order to apply the energy estimates of [15], we need that the coefficients of dthapo will 
be independent of t. This is because the specific form of the inner-product in Hg ^ spaces 
which takes into account the matrix of the system (3.1). In Section 5 we will further 
clarify this issue. Therefore we divide the second row by —g^^ and in order to preserve 
the symmetry of the system, we also multiply the third row by {—g^^)~^ . Thus the wave 
equations (2.5) are equivalent to the system 



(3.17) 



dtgafi 

dthal30 



h 



al30 



[-g^'^y {2g^^dah^^o + g'^'daK/sb + H^/s{g, dg) 
Sttw^ ((1 - Kw^)gai3 + 2(1 + Kw'^)UaUp) 



-g'T'g^'dtKpa = {-g'T'g'^'daKpo 

To shorten and simplify notations, we introducing the auxiliary variables 

(v, dtv, d^v) = {gap - Vap, dtgap, d^gap), 

where rjap denotes the Minkowski metric and dx denotes the set of all spatial derivatives. 
We also set (cq)" = (1,0,0,0) and W = { w, u°' — Cq) stands for the Makino and the fluid 
variables. Finally, 

U = {v, dtv, d^v, W) 

represents the unknowns of the coupled systems. 

We write the matrices in a block form, A = (ajj), the k x k identity matrix is denoted 
by and Omxn is the zero matrix. 

The coupled system (3.17) and (3.8) can be written in the form of (3.1), where A°' are 
55 X 55 symmetric matrices which depends only on v and W. We shall describe now the 
structure of these matrices: 



(3.18) 



where 



A%v,W) 



( eio 

10x10 OioxSO Oiox5 

OlOxlO Gio 040 x 30 Oiox5 

O3OXIO O30XIO ^^-^ 



y OsxlO UsxlO 



0, 



*33 
OsxSO 



30x5 
1-44 



*33 



lOO 



^^^eio ^^^eio ^^^eio 
^2ieio ^^^eio ^^^eio 
^^^eio ^^^eio ^^^eio 
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and = SLl^{ga/3,w,u°') is given by (3.8) when z/ = 0. From (3.17) we see that the 
coefficients of daU, a = 1, 2, 3, have the form 



OlOxlO 


0l0x40 


OsxS 


040X10 




O4OX5 


OsxlO 


05x40 


„a 
•^33 



where agg = a^^{ga/3,w,u°') is from the system (3.8) of the fluid and 

\ 
/ 

An essential demand is that 3.22(9013) G Hg s, whenever gap — rjap € Hg s- Obviously, this 
does not hold for the matrix in (3.19). Therefore we need to modify these matrices by a 
constant matrix 

/ OioxlO 









fl-^^eio g'^^eio g^^eio 


(3.19) 






030 x 30 



-22 



'We 



then {a.22 — C22) (v) G Hg s whenever v G Hg g. So we set 



<5«ieio 6-^e,o S^^e.o \ 







30x30 



(3.20) A''{v,W) 
and a constant matrix 



10x10 
40x10 
5x10 







10x40 







'■22 

5x40 



-22 



05x5 
O4OX5 







( OioxlO 


0l0x40 


05x5 


(3.21) 


r = 1 


040X10 


^22 


O4OX5 






\ OsxlO 


05x40 


05x5 



We turn now to the lower order terms. The presence of the fractional power w'^^^"' in 
(3.17) causes substantial technical difficulties. We set 



(3.22) 



f{v,W) 



((1 - Kw^)gap + 2(1 + Kw^)uaup) 



then we can write B{U) in the form 

B{U) = B{U){v, dtv, d^vf + F{v, W), 
where J^{v, W) = (0, f{v, W), 0, 0)^ and 



(3.23) 



BiU) 



^ OlOxlO 


eio 


OlOxlO 


OlOxlO 


OlOxlO 


\ 


OlOxlO 


b22 


b23 


b24 


b25 




030X10 


O3OXIO 


O3OXIO 


O3OXIO 


O3OXIO 




\ 05X10 


b42 


b43 


b44 


b45 


/ 
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The block h2j, j = 2, 3, 4, 5, appear from the quadratic terms in (2.4) 



where C^^Ti^f' are a combination of Kronecker deltas with integer coefficients. Thus 



The block h^j, j = 2,3,4,5, appear from the multiplication of the reflection F^^ and u'' 
in (3.8) with the Christoffel symbols. So its coefficients consist multiplications of Qa/s, g"^^ 
and . 

In summary, we can write the coupled systems (3.17) and (3.8) as a symmetric hyperbolic 
system 

(3.24) A\v, W)dtU = i{A%v, W) + dJJ + B{U) {dtv\+ T{v, W), 

\dxvj 

where A^{U) is positive definite in the neighborhood of the initial data (2.6), A^{0) —655 = 
^"(0) = and is a constant symmetric matrix. 

4. The Hs^s spaces and their properties 

The definition of the weighted Sobolev spaces of fractional order Hg^s, Definition 2.1, is 
due to Trieble [25]. Here we quote the propositions and properties which are needed for 
the proof of the main result. For their proofs see [4, 19, 25] 

We start with some notations. 

• Let {ipj} be the sequence of functions in Definition 2.1. For any positive 7 we set 

00 

j=0 

and we will use the convention HmHh^ 4- ^ = 4- The subscripts 2^ mean a scaling 

by 2^ that is, {^]u)2j{x) = (^/^Jm)(2%). 

• For a non-negative integer m and /3 e M, the space CJ^ is the set of all functions 
having continuous partial derivatives up to order m and such that the norm 

(4.2) \\u\\c^ = V sup ((1 + |x|)^+l°l|a°M(x)|) 

|a|<m, 

is finite. 

• We will use the notation A < 5 to denote an inequality A < CB where the positive 
constant C does not depend on the parameters in question. 

Proposition 4.1. For any positive 7, there are two positive constants 09(7) and £1(7) such 
that 



(4.3) 



co(7)||m||h,., < 



,(5,7 



< Ci(7)I|m||//, 



16 U. BRAUER AND L. KARP 

Proposition 4.2. For any nonnegative integer m, positive 7 and 5 there holds 

(4-4) II^IIL,^,. ^ < IkllL,.,.' 

where ||^^||m,5 is defined by (1-7). 
Proposition 4.3. If si < S2 and 61 < 62, then 

(4-5) h\\H,,.6, < M\h,,,s,- 

Proposition 4.4. If u G Hg^^, then 

(4-6) < \\u\\h,,,- 

Proposition 4.5. Let si, S2 > s, Si + S2 > s + |, si + S2 > and 5i + 82 > 5 — ^. If 

u G Hs^^Si and v G -f^sj.fc' ^^^^^ 

(4-7) ^ Mh... ,. \Hh.. 



Remark 4.6. If for a fixed constant Cq, m — Cq G and t> G -ffsa.fe) then we can apply 

the multiphcation property (4.7) to {u — Cq)v and obtain 

(4.8) < - Co||^^^^^^ + |co|) \\v\\h^^^,^ ■ 
Proposition 4.7. Let u G Hs^s fl l</3, 0<s</3 + i anc? 5 G M, i/ien 

(4.9) lll^l^b,,, <C(h|U^)||«||^,,,,. 
Proposition 4.8. Let s' < s and 6' < 6, then the embedding 

(4.10) Hs,s H,',5'. 
is compact. 

Proposition 4.9. If s > | + m and 5 + | > /3, then 

(4.11) \\u\\c^ < \\u\\h^^s, 
where \\u\\c^ is given by (4-2). 

Proposition 4.10. Let F : R™ — j- M} be C'^^^ -function such that F(0) = and where 
N > [s] + 1. Then there is a constant C such that for any u G Hg^s 

(4.12) WmWn^, < C\\F\\cN,. (1 + \\uU^,. 

Proposition 4.11. 

(a) The class C^{M?) is dense in Hs^s- 

(b) Given u G Hg s, s' > s > and 5' > 5. Then for p > there is Up G C^(M^) and a 
positive constant C{p) such that 

(4.13) \\up-u\\H,,<P and \\up\\H^,^,, < C{p)\\u\\h,,s- 
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4.1. Product spaces. The unknown of system (3.24) is a vector valued function 

U = {v, dtv, d^v, W) , 

where v = gap — rjap stands for the fields variables and W = {w, — Cg ) stands for the 
fluid variables. We consider it in the space 

(4.14) Xs^s := Hs^s X -f^s,<5+i x Hg^s+i x Hs+1,5+2, 
with the norm (see (4.1)) 

(4.15) \\U\\' = \\v\\' + WdMll ^^, ,, + ll^x^llk.... + II^IIL..... • 



Remark 4.12. Note that if f/ G Xg^s, then v G Hg^s and dxV G Hg^s+i, so by the integral 
representation of the norm Hs^s (see [1, §2)]), we have that 




5. Energy estimates 

In this section we will derive the energy estimates for a linear symmetric hyperbolic 
system which we have obtained by linearising (3.24). So we consider 

(5.1) A^dtU = {A"" + r)daU + B 

where U = {v, dtv, dxV, W), the matrices A*^, A"", B and have the same structural form 
as the corresponding matrices in (3.24), is a constant matrix, and the vectors J-" and V 
have the form (0, /, 0, 0) and (0, d2, d^, d^) respectively. 

Assumptions 5.1. All the matrices have the same block structure as (3.18), (3.20) and 
(3.23) and satisfy: 

(5.2a) {A%t, ■) - ess) , ^"(t, ■) e Hs+1,5; 

(5.2b) 3 Co > 1 such that co'^V^V < V'^ A^V < cqV^V, W G M^^ 

(5.2c) 9t^°(t,-) G 

(5.2d) h2j{t,-)eH,,s+u J = 2,3,4,5; 

(5.2e) b4,(t,-) e J = 2,3,4,5; 

(5.2f) T{t,-),V{t,-)eHs,5+i. 

5.1. Xs ,5-energy estimates. We turn now to the definition of an inner-product in the 
space Xg^s which takes into account the structure of the matrix A^ of the system (5.1). 
Let F{u) denote the Fourier transform of a distribution u and set 

A^(m) = (1 - A)i = ((1 + leP)^ F 

The standard inner-product on the Bessel-potential spaces is 

{Ui,U2)s = (A'(Mi), A'(m2))l2- 



U 



18 U. BRAUER AND L. KARP 

Taking into account the term- wise definition of tlie norm (2.1), we define the inner-product 
in s as follows: 

oo 

(5.3) («i,«2)s,5 := 5^2(^+1)2^- (a^ (V^l^i),, , {i^]u,) 



i=o 

where {u)2j denotes scaling by Hence {u,u)s^s = W'^lfn^g^ — W'^W'h^ g by Proposition 4.1. 
To each component of the space 

we assign its own inner-product. Since A'^ = (a°j), where a°- is the zero matrix for i ^ j, 
a° is the identity for i = 1,2, we assign to the first two components the inner-product 

(5.3) , while for the other terms we insert term-wise. 

Definition 5.2 (Inner-product in Xs.s). Let Ui = {vi,dtVi,dxVi,Wi) G Xg^s, « = 1,2 and 
assume that the matrix A^ satisfies Assumption 5.1, then we have 

• Inner-product on Hg s: {vi,V2)s,s, where is defined by (5.3); 

• Inner-product on Hg^s+i- {dtVi,dtV2)s,5+i, where is defined by (5.3); 

• Inner-product on Hg^s+i- 

oo 

(5.4) (5.^i,5.^2).,,+,,.o^ := 5^2(^+^+t)2^- (^|9.^02. ' (^33)2. 

j=0 

• Inner-product on 

00 

(5.5) (W^i,^2).,,,,,.o^ := j:2i^-^-i>^ , (a°,),, A-^ 

• Inner-product on Xg^s: 

(5 6) 

We denote by \\U\\y the norm associated with the inner-product (5.6). By assump- 
tion (5.2b) the following equivalence holds 

(5-7) \\U\\y <\\U\\y <\\U\\x ■ 

In order simplify the notation we set U{t) = U{t,x^,x'^,x^). 

Lemma 5.3. Let | < s, 5 > — | and assume the coefficients of (5.1) satisfy Assumptions 
5.1. IfU{t) e C^{R^) IS a solution of (5.1), then 

(5.8) I (f/(t), t/(t))x^,,,,o < Ceo ((f/(t), t/(t))x + 1) , 

where the constant C depends on the corresponding norms of the coefficients, s and 5. 



{Ul,U2)x , ,n := (^'1,^2)^,5+ {dtVl,dtV2)s,5^ 
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The corresponding energy estimates for the vacuum Einstein equations were obtained 
in [15]. The same techniques can be apphed here with some obvious modifications. We 
therefore give only a short sketch of the proof. 



21 



L2 



Sketch of the proof. From the inner-product (5.6) we see that 
1 d 

oo 

oo 
j=0 

By the Cauchy Schwarz inequahty, we obtain 

and by Assumption (5.2c) the first infinite sum is less than 



2J 



L2 
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3=0 

II „.l|2^ 

's,<5+l,2 



L2 



=^ll^*^33Loo \\dxV\\H 

A similar estimate holds for the second infinite sum 
The most difficult part is the estimate 



(5.9) {dtv, v)s,5 + {dxV, d^dtv) 



,5+1, a' 







^ \\u\\\ 



and here it is essential to use the assumption that G Hs+i^s and s > |. We refer 
to [15, §4] for the details of this estimate. The estimate of {W, dfW) a'^ relies on 
similar ideas to those of (5.9), but it is simpler, since the regularity of the corresponding 
coefficients of (5.1) and is (s + 1). 

Having collected the estimates of all the terms, we have 



-^(U,U)y <\\u\\l 

and by the equivalence (5.7) we obtain (5.8). □ 

Remark 5.4. In order to use the fact that U{t) satisfies system (5.1) in the proof of Lemma 
5.3, we have to replace the expression aggA* (^tjj'jdtdxV^ which appears in the inner-product 
(5.4), by A^ (^ip'^ a'^^dtdxvY This is achieved by using the Kato-Ponce commutator for [24, 
§3.6] in several steps. In one of them we apply it to the commutator 

(5.10) i^A),. (A^5.) i^ptv),, - (A'd^) i^Ai^ptv), 



'2i 
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where = ipk (^ji^j^ ■ Since A'^dx is a pseudo differential operator of order s + 1, 
it follows from [21, §3.6] that the norm of (5.10) is bounded by (^ipjdtv^, 

ll(*fc(a!]3-e3o))2.L., see [15, §4]. 

However, had the coefficients ofdt{dtv) be depended on t, the inner-product {dtVi, dtV2)s,s+i 
in (5.6) would have been replaced by 



'2J 



and 

H" 



3=0 



21 



L2 



for a certain 10 x 10 matrix a^2- And now we would have needed to make the same 
commutation as above, but with the operator h!^dt instead oi K^dx- Since h.^dt is a pseudo- 
differential operator of order s, this required that s — 1 > | and we would not have gotten 
the desired regularity. We therefore transfered the system (3.16) to (3.17). 

Note that for the inner-product (5.4) we use the operator for dtv, dxV, while for the 
fluid variables W, we have the operator A*+^ in the inner-product (5.5). Therefore, for the 
Euler equations we can perform the operation mentioned above under the condition s > |. 

5.2. L^-energy estimates. The L^(M'^)-space is the closure of all continuous functions 
under the norm 

II^IIl^ = / (1 + 1^1)^ \u{x)\'^dx. 
This norm is equivalent to the norm ||tt||_H-o^ (see [25]). Similar to (4.14), we set 

and \\U\\l = Ml. + \\dtv\\l.^^ + \\dxv\\l.^^ + \\W\\l.^^. 

We also define the inner- product on Ys in accordance with to the system (5.1): 



{UU U2)y„A^ " y + \x\f'v'[v2dx + y (1 + k|)'('+'^ (dtV.f {dtV2) dx 

+ J{l + \x\y^'+'^ {dxvif43 {d.V2)dx + j{l + \x\y^'+'^W^aP,^W2dx 
and the associated norm ||t^||y^_4o = {U,U)ys,ao- 

Lemma 5.5. Assume the coefficients of (5.1) satisfy Assumptions 5.1. IfU{t) E Xi^s is a 
solution of (5.1), then 



d_ 
It 

where the constant C depends upon the L°° -norms of A", daA" and B. 



- {U{t),U{t))y^^^. < CC, [{Uit),Uit))y^^^, + + \\V\\l.^^ 
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6. The iteration process 

In this section we adopt Majda's iterative scheme [17] in order to prove the well-posdeness 
of the coupled hyperbohc system (3.24) in the if^^^-spaces. A similar approach was carried 
out in [15] for the vacuum Einstein equations, but in the presence of a perfect fluid there 
are additional difficulties. Since the density is not strictly positive, the zero order term B 
which contains the fractional power of the Makino variable ly^/^T"^) is not necessarily C°°. 
Hence we could not apply the standard existence theory for symmetric hyperbolic systems. 

We denote the initial data by Uq = (0, ip, dx(f>, Wq), where Wo = {wq, {uq)°' — Cq) repre- 
sents the initial data for the Makino variable w and the four velocity vector m° — Cq. 

Theorem 6.1. Let | < s < + > ~|- Assume Uq G Xg^s, Wq > and that there 
exits a positive constant fi such that 

(6.1) -V^V < V^A\(p, Wo)V < fiV^V for all V G M^^ 

Then there exists a positive constant T and a unique solution U{t) = (f (t), dtv(t), d^vit), W{t)) 
to the system (3.24) such that U{0,x) = Uo{x), 

(6.2) U eCi[0,T],Xs,s) and W E C {[0,T], Hs,5+3) ■ 

The proof of Theorem 6.1 will be carried out in several steps: 

1. Setting up the iterative scheme; 

2. Proving that the fractional power (t(;'=)2/(7-i) jg ^ C^-function; 

3. Boundedness of the iteration sequence in the X^^^-norm; 

4. Weak converges of the iteration sequence in the Xs^s norm; 

5. Uniqueness of the solution obtained; 

6. Continuity in the corresponding norm. 

A part of the above proofs are standard, but some of them require a special attention 
due to the specific form of the system (3.24) and use of the space X^ s. Moreover, the fact 
that the matrices A'^ = ^"(f , W) are not dependent on the derivative of the semi- metric 
plays an essential role here. 

Step 1. From condition (6.1) and the embedding into the continuous. Proposition 4.9, we 
see that there is a bounded domain G C M^^ containing the initial value Uq and a constant 
Co > 1 such that 

(6.3) —V^V < V^A^h, WW < cqV^V 

Co 

for all U = {v, dfV, d^v, W) G G and V G M^^. By means of the density properties of Hg s, 
Proposition 4.11, there is a sequence 

(6.4) {u^}Zo = {(0^/,9.0^4,K)' -eo)}r=o ^ 

and a positive constant R such that 

(6-5) P'o\U,.,,<C\\Uoh^^^, 
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(6.6) \\U-U^\\y <R^UeG 

and 

DO — 

(6-7) l|Co'-CoL.,,<-^. 

The iterative scheme is defined as follows: let U^{t,x) = Uq{x) and U'^~^^{t,x) = 
{v^~^^{t,x),dtv''~^^{t,x),dxv''~^^{t,x),W'^~^^{t,x)) be a solution of the linear Cauchy prob- 



lem 



(6.8) 



A'^{v'',W'')dtU''+^ = {A^{v'',W'')+C'')daU''+^ 



f/fc+i(0,x) = U^+\x) 

where J^(w^ ly'^) = (0, /(w^ VT'^), 0, 0), /(^;^ W'') is given by (3.22) and W'' = (M;^ (u")^- 



Step 2. The iterative method relies on the fact that solutions of linear symmetric hyper- 
bolic systems with coefficients and initial data, are also C^. However, even if ■u;*' > 
and E C^, it does not guarantee that (u7'^)2/(7-i) jg ^ C^-function. Since the function 
/(f'"', W'^) contains the term (w''y^^"'~^\ we must assure that it is a C^-function. 

Proposition 6.2. Let u G Hs^s he non-negative and /3 > 0. Then there is a sequence 



{n^'} C such that u in the Hs s-norm and {u^Y ^ C*, 



oo 

• 



Proof. Let £ > 0, then by Proposition 4.11 there is G with \\u — u'^\\h s ^' '^^^^ 
now a positive number M so that supp(M^) C < M}, and let Xm be the cut-off function 
satisfying Xm{x) = 1 for \x\ < M and Xm{x) = for |a;| > M + l. For any positive number 
g we set 

u'^'^ix) = xnix) {u%x) + g) . 

Then {u'^'^Y ^ C*^, since {u^ + g) > Q and xm is a cut-off function. Moreover, u^'^ — = 
XmQi hence u^'^ — > in the if^.^-norm as — 0. □ 

Thus we may assume that {wolfc^c approximation of the initial data of the 

Makino variable wq in (6.4), satisfies (m;^)2/(t-i) g Co°°. We turn now showing that for 
t > 



e\t,x) = {w'')--^ (t,x) 
Proposition 6.3. For each integer k >0, e^'(t, ■) G C^(M'^) 



is also a C^-function 
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Proof. We conduct the proof by induction. Obviously the statement holds when k = 0. 
Now the SI*'* row of the system (6.8) is equivalent to (3.5a), so the linearization of it results 
in 

where the -?„(■, ■) is the projection in equation (3.5a) and (r^^) are the Christoffel symbols 
with respect to the semi-metric (7^^. Since by [1, Theorem 2.6], u^^O, x) > 1, we can divide 
(6.9) by (u^)^ and we conclude that e'^'^^ satisfies a first order linear equation of the form 

5te^+i + ba{t, x)dae^+^ + c(t, x) = 
e'=+i(0,x) = (^/^^^^(x))^ 

Note that c(t,x) contains the term e^, but this term is by the induction hypothesis. 
Hence all the coefficients of (6.10) are C^-functions. 

We solve (6.10) by means of the characteristic method. So let $(s,y) be the solution of 
the system 

^ = 1, ^ = 6,(t,x), t(0) = 0, x(0)=y. 
as as 

Then e*'"'"^(t, x) = Z($~^(t, x)), where Z{s,y) is the solution of then initial value problem 

rl7 2 
— = -c(.,x), Z{0,y) = {w',^\y)) — . 

Obviously 

Z{s,y) = {w^+\y))^ - I c{T,x{r))dT. 



(6.10) 



Since c(r, ■) G and (w;o+^)^-i G and by Proposition 6.2, we have that Z{s, ■) e C^ 
and hence also (w^'+i)^ (t, ■) = e'^+^t, ■) G C;^. 



00 




□ 



Step 3. We conclude from Proposition 6.2, and the theory of linear symmetric hyperbolic 
systems (c.f. [11]), that for each k there is a solution U^{t,x) of the linear system (6.8) 
such that U''{t, ■) G ^^(M^). Therefore by (6.5), (6.6) and (6.7), for each k we have 

Tfc = supjr: sup ||f/^'(^) - f^oL <^1>0. 

t 0<t<T J 

Proposition 6.4. There are positive constants T* and L such that 

(6.11) sup It : sup \\U\t)-U^\\^ <r\>T* for all k 

[ 0<t<T '''^ J 

and 

(6.12) sup ||5tW^*^|L <L for all k. 

Q<t<T* ^'iis,S+3 
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Proof. We prove it by induction. Set = [7^+^ — U^, then 
(6.13) 



where 

(6.14) = [A^iv'', W'') + r) daU^ + i3(f/^' 





and V^+^{{),x) = U^^\x) - U^{x), 

In order to apply the energy estimate, Lemma 5.3, to (6.13) we have to check that 
Assumptions 5.1 are satisfied and the corresponding norms are independent of k. Clearly 
all the matrices have the same structure. From (6.3) we see that condition 5.2b holds. By 
the induction hypothesis, 

h'(t) - + \\d.v\t) - d.Al\f^^^^^^ < \\U\t) - u^,\\l^^ < R\ 

therefore by Remark 4.12, |k'^'(t) — ^oHh- R- Applying the Moser type estimate. 

Proposition 4.10, we have 



\A\v\ W^) - 65511' ^ < { V\t . + IIW^'"' 



sl|c'-Co°iiL + li^o°iiL^(«^+iwiiL 



In a similar way we get that ||^"(f VT'^) ||^ is bounded by a constant depending on 
R. By Propositions 4.3, 4.5 ,4.10, Remark 4.6 and the structure of the matrix B in (3.23), 
we obtain that ||b2j(f/^) ||^ < ^ (-R + l|f^oll^si) ^^"^ ^ similar estimate holds 

for ||b4,(f/'=)||^^^^^^, j = 2,m!5. 

We recall ihe J^{v^ ,W^) = (0, /(^;^ ly'^), 0, 0), where f{v'',W'') is given by (3.22). Ap- 
plying again Propositions 4.5, 4.10 and Remark 4.6, we obtain that 

ll/{"''.»"')L.,„. ^ Ih'f^lL (ll-i^M,. + IK"")' -'Sll«.,„. + 1) ■ 

Now, for I < s < + |, we apply the estimate for the fractional power. Proposition 4.7, 
and together with Propositions 4.3 and 4.9, we have 



Hs,s+i 

Using similar arguments as in the previous estimates, we conclude that ||/('y'^, 

is bounded by a constant depending on R but not on k. 

The required estimate for follows from the multiplicity property (4.7) and the esti 
mates which we already obtained for L4"(f*^,iy'^) „ and LSff/'^) „ 
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It remains to verify (5.2c); note that by the induction hypothesis (6.11), condition (6.6) 
and the embedding (4.11), we have 

9^° 



sup 

G 



\dtA%v^,W^ 



< 



sup 

G 



dv 



dW 



dtW^\ 



< 



:<5+l 



Since 11 ^VT^ 



is also 



is bounded by hypothesis (6.12), we see that 9t^°(w'^, W^) 

bounded by a constant depending on R but not on k. 

We can now apply Lemma 5.3, and with the combination of Gronwall's inequality, con- 
dition (6.7) and the equivalence (5.7), we conclude that there is a constant C depending 
on R such that 

2 ^ ^CcoT I R 



sup \\V^{t)\ 

0<t<T 



< e 



CclT 



Hence T* = sup{T : e^'^"^ (^f- + Cc^t) < R^} > 0. 

We show now (6.12). It follows from the structure of the matrices A^, A"" and B (see 
(3.18), (3.20) and (3.23)) that 

(6.15) dtW'+' = {4,{v\ W'))-' [4,{v\ W')daW'+' + h,^,^,){U')dy+'] . 



Since the Wv'^ 



+ 1,S1 



and ||?^'^(t)||xs 5 are bounded by a constant in- 
dependent of k, it follows from (6.15) and by means of calculus in the if^^^-spaces that 



s+l,i5+2 



fc+il 



,5+3 



is also bounded by a constant independent of k. 



□ 



Step 4. Having shown the boundedness of U'^ in 5-norm, we use the compact embed- 
ding. Proposition 4.8, and conclude that there is f/ G Xg'^s' such that ||?7'^(t) — ?7(t) 11 — 

for any s' < s and 6' < 6. Furthermore, by Remark 4.6, v''{t) — ?■ v{t) in if^z+i ^/-norm. 
Thus if we choose | < s' < s and — | < 5' < 5, then by the embedding (4.11), 

v''{t) v{t), W''{t) W{t) in 

and 

dtv''{t) dtv{t), dtW'it) dtW{t) in 
Thus, U{t) = {v{t),dtv{t),d:,v{t),W{t)) is a solution of the system (3.24). 

Proposition 6.5. For any $ G Xg^s, 

(6.16) l\m{U'it),^)^^^ = {Uit),^)^^^ 

uniformly for <t< T* , and where (■, ■)-^ ^ denote the inner-product (5.6) with JP being 
the identity matrix. 
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As a consequence of the weak convergence (6.16), we have 

l|f^Wllx.,<limmf||f/^t)L^,. 
For the proof of Proposition 6.5 see [15, §5]. 

Step 5. Here we shall prove uniqueness. Our method relies on the L^-energy estimates, 
Lemma 5.5. A similar approach was used in [15] for the vacuum Einstein equation, we 
shall therefore emphasize the estimates of the terms which do not appear in the vacuum 
equations. 

Proposition 6.6. Suppose Ui(t),U2(t) G Xg^s are two solutions of (3.24) the same 
initial data, then Ui{t) = U2(t). 

Proof. Let V{t) = Ui{t) — U2{t), then it satisfies the linear system 

(6.17) A\vi, Wi)dtV = A^ivi, W,)daV + BiUi)V + ^(t;i, W^) - T{v2, W2) + V, 

where 

V = - {A\v^, W,) - A\v2, W2)) dtU2 

+ {A\v^, W^) - A\V2, W2)) daU2 + (i3(f/i) - B{U2)) U2. 
Applying Lemma 5.5 to system (6.17), we have 
d 



^1+1 



Thus, our main task is to bound the two last terms on the right hand side of (6.19) by 

\\ui{t) - U2ml,- 

Staring with the first one, we have that 



11^(^1, W{) - Hv2. W2) = 11/(^1, W,) - fiv2, W2) 11^. 

0+1 0+1 



and 



f{v,,W,)- f{v2,W2) 

^ {tv^ + (1 - r) V2, tW, + (1 - r) 1^2)) (v^ - V2) 
^ {rv, + (1 - r) V2, tW^ + (1 - r) H^a)) (W^i - W2) 

for some r G [0, 1]. Note that vi,V2 belong to Lj, while we have to take the L^_,„-|^-norm 
of the above expressions. However, since the Makino w G Hs+i,s+2 C Hs,s+2, we get 
from (3.22), and Propositions 4.7 and 4.5 that |^ G Hs s+2- Therefore, by the equivalence 



LOCAL EXISTENCE OF RELATIVISTIC PERFECT FLUIDS 



27 



\u\\i2 ~ ||^t||//oa (Proposition 4.2) and the multiplicity property (4.7), we obtain 



< 



dv 
df 



^s+i 



dv 



H, 



0,(5+1 



dv 



{v1-V2)\\h,,, < 



S+2 



Of 



dv 



H 



[Vi-V2)\\l2 



The other term is somewhat easier to treat, since ||^||^oo is bounded by a constant 
depending on ||f/i||Y and ||t/2|lx ' hence. 



df 



dW 



{Wi - W2 



< 



\\{W^-W2)\\ 



s+2 



Thus 



\\f{VuW,)-f{v2,W2)\\l2 <\\U^-U2 



S+l 



In- 



We shall now estimate the first term of V in (6.18). From the structure of A^{v, W) we 
see that 

{A'{vi,Wi)-A'{v2,W2))dtU2 = 

{4,iv,) - 4,{v2)) dtd,V2 + {s-Uv,, W,) - aP,,{v2, W2)) dtW2, 
and furthermore, 



dal 



■33 



dv 



(tVi + (1 - t)v2) (Vi -V2) 



for some r G [0, 1]. At this stage we apply the calculus in the if^^^-spaces. The multiplica- 
tion property (4.7) is used twice, once with si = 1, S2 = s — 1 and s = 0, and the second 
time with si = s, S2 = 1 and s = 1. All these result in 



i^si'^i) - ^si'i^^)) dtd.^V2\ 



T2 



< 



< 



9a! 



33 



dv 



{rvi + (1 - t)v2) {vi - V2) 



(a°3(t;i) - a°3(t;2)) dtd^^V2\\ ^^^^^ 
\\dtdxV2\\ TT 



{tVi + (1 - t)v2) 

dv 



11(^1- ^2) iIh,,, ii^*^2|i^,_, 



< \\{rVi + (1 - t)v2)\\\ \\{Vi - V2)fH^JdtV2\\\^ . 



By similar calculations, 

UlA^l^Wr) - S.I^{V2,W2)) dtW2\ 



[2 ^ 



ds.^4 



dv 



X \\dtW2\\H, 



s.S+3 



[Vi-V2)\\h^^^ + 



dW 



11(1^1-1^2)11 



^S+2 
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We recall that ||f^tM^2|lH^ ^+3 is bounded by (6.12) and — t'2)||^^_^ ^ IK^i ~ ^2)11^2 + 
\\{d^vi - d^V2)\\l2 , hence, 

(5+1 

II W,) - A\V2, W2)) dtU2\\l,^^ <\\Ul- f/2||y, . 

The remaining terms can be estimated by similar methods. Thus we conclude from 
(6.19) and the equivalence \\V\\y ^ IIFIL that 

^l|fWllL.5l|vWllL.. 

and since V{0) = 0, Gronwall's inequality implies that V{t) = 0. □ 



Step 6. Since Xg^s is a Hilbert space it suffices to show that 
(6.20) limsup||f/(t)||^ < WUmix ,,0 

^_^Q+ s,d,A^' s,d,A^ 

in order to establish the continuity in the norm. Here depends on the initial data 
and Wq, that is, = A'^{4>,Wo). The proof of (6.20) relies on the same arguments as in 
[17] and we therefore leave it out. This completes the proof of Theorem 6.1. 

7. Proof of the main result 

The proof of the main result. Theorem 2.3, actually follows from Theorem 6.1, we just 
have to check whether the initial data of the gravitational fields and of the fluid satisfy the 
assumptions of Theorem 6.1. 

We recall that v{t) = Qapif) — rjap, so setting (p = v{0), we have by the assumptions of 
Theorem 2.3 that (j) G Hs+i,s, and by (2.6), ^.^^{(p) is a positive deflnite matrix. The initial 
data for dtv is given hj ip = dtQaisiO), dtQabi^) = -'^Kab {a,b = 1,2,3), where Kab is the 
second fundamental form and by the assumptions of Theorem 2.3 it belongs to Hs,s+i- The 
remaining data, dtQaoiO), are determined by the harmonic gauge (2.3). Since we assume 
that ((70/3(0) — rjap.Kab) G Hg+i^s x Hg^s+i, it follows from Propositions 4.5 and 4.10 that 
also dtQaoiO) e Hs,s+i- 

The construction of initial data for the fluid variable was studied in details in [4, Theorem 
2.6], from which it follows that Wo{-) = W{0,-) G Hs+i^s+2- Since is a unite timelike 
vector, it follows from Theorem 3.3 that a44(lVo) is a positive deflnite matrix. Therefore 
A^{(f),Wo) satisfles condition (6.1), and hence all the assumptions of Theorem 6.1 are 
fulflUed. 

So as a consequence of (6.2) and Remark 4.12, v{t) G C {[0,T], Hg+i^s). Hence (2.7) 
holds. Also, by (6.2), (2.8) holds too. That completes the proof. 

8. Appendix 

The classical paper of Hughes, Kato and Marsden [13] established the short time exis- 
tence of the vacuum Einstein equations by solving a second order quasi-linear hyperbolic 
system whose solutions {gap^dtQap) belong to W^^^ x H" for s > |. On the other hand. 
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Fisher and Marsden treated the Einstein vacuum equation by means of the theory of sym- 
metric hyperbohc system. However, they only obtained the regularity of s > |. 

In [15] we generalized the result of [13] to the Hg^s spaces, treating however, the Einstein 
equations as a symmetric hyperbolic system. Since the techniques of [15], and in particular 
the energy estimates, play an essential role in the present paper, we outline its main idea 
that enables us to obtain the same regularity as in [13]. 

We present a heuristic argument explaining the essential idea. First, if a function v 
satisfies a wave equation, then the vector V = {v, dfV, d^v) satisfies a symmetric hyperbolic 
system. The general condition for existence and uniqueness in the if^(M'^) spaces is s > |. 
Hence, we have by this method that dtv, d^v G H'^ for s > |. 

However, in our case we improve this regularity to {v, dtv, dxv) G H^'^^ x H'^ x for 
s > |. This is because we do not consider a general quasi-linear symmetric hyperbolic 
system where the matrices A'^iV) depend on V , but a system in which the matrices A'^{v) 
only depend on v but not on its derivatives. 

To see how this fact allows us to improve the regularity of the solution we will derive 
energy estimates for the linearized symmetric hyperbolic system. For the sake of clarity 
we consider a simple hyperbolic system 

(8.1) dtV = A-{v)daV, 
then its linearized form is 

(8.2) dtV = k'daV. 

Note that in each iteration we solve the linear system (8.2) with A"- = A"(f^), and since 
V'' = (w^ a^w^') G H% G H^^^ and hence A" = I"(w'=) G H'+^ by Moser type 
estimates. 

The crucial step is to derive the energy estimate 

(8.3) |(l||^f„.)<Cri|i. 

for s > I and whenever V satisfies the linear system (8.2). We recall that ||V^||_f/'i = 
||A*y||j;^2, where A'^ is the pseudodifferential operator (1 — A)t. 

One of the basic tools for obtaining (8.3) are commutator's estimates. Here we shall use 
the following Pseudodifferential operators version of Kato-Ponce estimate [21, §3.6]: Let 
P be a differential operator in the class OPS^ q, then 

(8.4) \\Pifg) - fP{9)h. < C{||V/|Uo.||^||^.-. + ||/||^,.||^?|U°=} , 
for any f e H' H and g e H'~^ n L°°. 
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The standard way to obtain (8.3) is to differentiate ||V^||^3 with respect to time, to insert 



the differential equation (8.2) and then apply a suitable commutator which leads to 

'2 {A'{V), A' (dtV))^, = { A'iV), A' { A^daV 
A^(\/),>(A^ (daV)) 



--iin 



L2 



+ (A^(V), A'(A''daV)-A^A'{daV)) 



L2 



and then the first term is taking care by integration by parts and the second one is by 
applying the above Kato-Ponce estimate to the operator A'^. 

But this procedure results in a term of the form ||9aV^||j;^oo which contains ll^aS^^f ||loo. 
In order to estimate it by ||i9a9xt'||i^s-i < ll^^f we need to require that s — 1 > |, and 
hence we do not get the desired result. 

We circumvent this difficulty by writing 

^daA'')v, 



(8.6) A'^daV = da(A''V 

and making the commutation 



A' (dj A^V 



{A' da) (A-v)- A'^ {A' da) {V) + A'^ (A'da) {V) . 



(8.7) 

which we insert into the first row of equation (8.5). Then we have to estimate three terms: 
(8.8) / = (a'{V), \{A'da) (a'^v) - A'' (A'da) {V) 



L2 



(8.9) 

and 
(8.10) 



II 



III 



A%V),A'' iA'da)iV) 



A%V),A'(idaA'')V 



L2 



L2 



For the first term we apply the Kato-Ponce commutator (8.4). However, this time we 
do it for the operator [A'^da) which has order s + 1, and hence 



(8.11) 



|/| < \\V\ 



< \\v\ 



(A'^da) [A-VJ - A'^ (A'^da) (V) 
\\VA^lo.\\V\\h^ 



L2 



So by Sobolev embedding theorem, we see that |/| < || || Likewise, since H'^ 

3 



is an algebra for s > ^, 



(8.12) 



< 



< WA"] 



Hs+i 



\V\ 



Since A^da = daA^ and A" is symmetric, we obtain a similar estimate for // by using 
integration by parts. Hence we conclude that the energy estimate (8.3) holds. Note that 
in the estimate of all three terms above we have used the fact that A"- G H'^^^. 
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For the general case where ^ I, one has to define an appropriated inner-product 
which takes into account the matrix A^. Details for the vacuum equations in the weighted 
spaces Hs^s and a positive definite A° can be found in [15, §4] and only slight modifications 
are needed in order to extend the energy estimates of [15] to the coupled system (3.24). 
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